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ABSTRACT 

We present and discuss the results of an extensive observational campaign devoted to 
GRB 071010A, a long-duration gamma-ray burst detected by the Swift satellite. This 
event was followed for almost a month in the optical/near-infrared (NIR) with various 
telescopes starting from about 2min after the high-energy event. Swift-XRT observa- 
tions started only later at about 0.4 d. The light-curve evolution allows us to single out 
an initial rising phase with a maximum at about 7 min, possibly the afterglow onset in 
the context of the standard fireball model, which is then followed by a smooth decay 
interrupted by a sharp rebrightening at about 0.6 d. The rebrightening was visible in 
both the optical/NIR and X-rays and can be interpreted as an episode of discrete 
energy injection, although various alternatives are possible. A steepening of the af- 
terglow Ught curve is recorded at about Id. The entire evolution of the optical/NIR 
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1 INTRODUCTION 



The Swift satelhte ( Gehrels et aL 2004 1 has revolutionised 



gamma-ray burst (GRB) science by introducing dedicated 
observational strategies and providing a wealth of high- 
quality data. In particular, a major breakthrough has been 
the spectral and temporal coverage at X-ray energies from 
a few minutes up to weeks or months after the GRB. This 
has provided unprecedented information on the late prompt 

Opti- 
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and early afterglow evolution (Nousek et al. 
cal observations from the ground, however, have not always 
been of comparable quality and late-time phases have been 
seldom monitored accurately ( Covino et al.|20^ |. This is a 
consequence of the strong increase in the number of events 
to be followed, the large average redshift of Swift GRBs 
( Berger et al.|2006 Jakobsson et al.|2006 1 , coupled with the 
difflculties in performing productive follow-up observations 
from the ground. While the early-time afterglow is often well 
sampled by the network of small and medium-class robotic 
telescopes around the world, the late-time afterglow can only 
be monitored with the largest instruments and only at the 
cost of significant observing time (Dai et al. 20071. Based 
on these considerations, we started a new program of af- 
terglow observations with the explicit goal of providing full 
coverage of the optical/near-infrared (NIR) afterglow light 
curve for a selected number of events (building on visibil- 
ity and availability of early-time observations) rather than 
providing sparse datapoints for most of the observable Swift 
GRBs. 

The importance of dense temporal and spectral sam- 
pling is particularly crucial for the identification and inter- 
pretation of the different features in the light-curve evolu- 
tion that are observed both in the optical (|Kann et al.|2007| ) 
and at higher energies ( Zhang fc Meszaros|2004 1 . In partic- 
ular, in spite of intense observational efforts and continu- 



ous theoretical interest ( 


Burrows & Racusin 


2006| 


Covino 


et al.| 


20061 ICurran et al.||2007| |Dado et al.||2007| |Dai et al. 


2007 


Panaitcscu 2008] e.g.,), the search for light-curve jet 



breaks fully satisfying the requirements of the standard fire- 
ball model, and thus allowing a correct inference of the true 



energy content of a GRB, is still far from complete (Liang 
|et aL|2007| |. 

The GRB071010A afterglow evolution was followed in 
the optical by TAROT and in the NIR by REM and GROND 
at early times, then by Gemini-N, GROND, NOT, TNG, 
NTT, Keck-I, Sampurnanand, and VLT in the subsequent 
two weeks. Multi-band observations were acquired to study 
the temporal evolution of the afterglow spectral energy dis- 
tribution (SED). A polarimetric observation was carried out 
almost in coincidence with a slope transition in the light 
curve. Swift XRT and BAT data were also studied. In Sect .[2] 
we report information about GRB071010A, in Sect.[3]our 
observations, reductions, and data analysis are described, 
and in Sect.|4]a full discussion is presented. Our main con- 
clusions are summarised in Sect.[5l 



2 GRB 071010A 
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GRB 071 010A was detected by Swift on Oct. 10, 2007 at 
03:41:12 ( [Moretti et "aLl|2007| ) (UT dates are used through- 
out this paper). The BAT light curve was relatively broad. 
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with a Tgo duration of 6 ± 1 s (a moderate-duration long 
GRB). The fluence was ( 2.0 ± 0.4) x 10~^er gcm~^ (un- 
certainty 90% confidence, Krimm et al. 2007). Automatic 
slewing to the GRB position with the Swift narrow-field 
instruments was disabled, so the XRT observations were de- 
layed by 34 ks ( Guidorzi et al.|2007b |. The optical afterglow 
was promptly identifie d by TAROT (qj2oo o = 19:12:14.624, 
^J2000 = -32:24:07.16; |Bloom et al.|2007[ ), revealing an ini- 
tial flux increase up to about 470 s after the high-energy 
event, followed by a regular decay (Klotz et al. 2007a|b I 
Spectroscopy was carried out with Keck (Prochaska et al 



2007), revealing a strong Mgll absorber and corresponding 
Fell lines at redshift z ~ 0.98. A refined analysis was car- 
ried out by|D'Elia et al.| (|2008[), deriving z = 0.985 ± 0.005 



(Icr error) for the absorbing system with the highest red- 
shift. Lines of Mgll (2796-2803 A), Mgl (2852 A), Sill 
(2336 A), Fell (2367, 2374, 2586, 2600 A), and Mnll (2594, 
2606 A) were identified. The afterglow was not detected in 
the radio band at 8.46 GHz almost two days after the GRB 
( Chandra fc Frail]|2007 l. Late-time observations led Perley 
et al. ( 2007 1 to suggest that the afterglow decay underwent 



a (possibly achromatic) steepening about one day after the 
burst. Galactic dust absorption along the GRB line of sight 
is Eb-v « 0.098 mag ([Schlegel et aL]p98l). 



Throughout this paper, the time decay and energy spec- 
tral indices a and /3 are defined by F{t, v) on {t — to)~°'i'~^ , 
where to is the trigger time of the burst. We assume a ACDM 
cosmology with = 0.27, = 0.73, and ho = 0.71. At 
the redshift of the GRB the luminosity distance is 6.6 Gpc 
(~ 2.0 X 10'^* cm), corresponding to a distance modulus 
H — 44.1 mag. All uncertainties are la unless stated oth- 
erwise. 



3 OBSERVATIONS, REDUCTIONS, AND 
ANALYSIS 

3.1 Observations and reduction 

The afterg low of GRB 071010A was observed by the TAROT 
telescope ( Bringer et al. 1999 1 with the white and red fil- 



ters, by the REM telescope ( Zerbi et al.||2001|l equipped 



with the REMIR NIR camera (Vitah et al. 2003 Conconi 



with GROND | |Greiner et "aL][2008| starting from 124 s af- 
ter the GRB (86s after reception of the GGN alert). The 
observations began at an initial airmass of about 2.3 with 
the target already setting. In order to obtain an accept- 
able signal-to-noise ratio (S/N) we needed to bin the REM 
data, thus decreasing the highest time resolution afforded 
by the instrument. Later NIR observations were carried out 
with Gemini-N equipped with NIRI, the TNG equipped with 
NIGS, and the NTT equipped with Sofl. In the optical the 
afterglow was again observed a few hours after the GRB 
with the Keck-I telescope equipped with LRIS, the Sam- 
purnanand telescope, the NOT equipped with PolCor, and 
the VLT with FORSl and FORS2. Polarimetry was carried 
out with the VLT equipped with FORSl on Oct. 11, 2007 
at 1:11:33 (0.89607d after the GRB) with a total exposure 
time of 1800 s. 

Data reduction was carried out following standard 
recipes with the Eclipse package ( Devillard|1997 1. Photom- 



etry was computed by aperture and profile-fitting photome- 
try with the Daophot ( Stetson] 1987 1 and SExtractor ( Bertin 



fc Arnouts|1995 1 packages. Photometric calibration was de- 



rived from the observation of standard stars in the optical on 



two different nights and with the 2MASS catalogue (Strut- 



skie et al. 2006 1 in the NIR, although in a few cases the 



small field of view of NIR detectors made the calibration 
more difficult. TAROT data were calibrated with observa- 
tions of two field stars and final magnitudes are expressed 



in the R2 system of the USNO-Bl catalogue (Monet et al 
2003 1 . The results from the photometric observations are 



reported in Table[T] 

For polarimetry, the offset between instrumental and 
intrinsic polarization was fixed by the observation of a po- 
larized standard star. Instrumental polarization was derived 
by observation of non-polarized stars. The procedure we fol- 
lowed for the analysis of the FORSl polarization data was 



extensively discussed by Covino et al. (1999 2002 20031 



At the time of the observation the afterglow optical emis- 
sion was consistent with null polarization with a 95% upper 
limit of 0.9% (1.3% at 3 a) for the degree of polarized flux. 

The Swift XRT data were reduced using the xrtpipeline 
task (v. 0.9. 9), applying standard calibration and filtering 
criteria, i.e., we cut out temporal intervals in which the CGD 
temperature was above — 47°C and removed hot and flick- 
ering pixels. An on-board event threshold of ~ 0.2 keV was 
applied to the central pixel; this was proven to reduce most 
of the background due to the bright Earth and/or CCD 
dark current. We selected XRT grades 0-12 for data all in 
PC mode. A circular region of 20 pixel radius was selected 
for source extraction. Physical ancillary response files were 
generated using the task xrtmkarf. Optical/NIR and X-ray 
light curves are shown in Fig. [2] 

The Swift BAT data were also analyzed. The 15- 
150 keV mask- weighted energy spectrum was extracted with 
the ftool hathmevt using the BAT refined position ( |Krimm| 



et al. 20071. All the required corrections were applied: we 



accounted for the slewing through the ftool batupdatephakw 
and produced the detector-response matrix with batdrmgen. 
The spectrum was corrected for systematics depending on 
energy with the ftool batphasyserr and was finally grouped 
by imposing a 3a threshold on each grouped energy channel. 



et al.|2004| , and by the 2.2 m MPI/ESO telescope equipped ^ ^ Analysis 



The optical/NIR light curve is complex, with an initial 
rising phase, a maximum, then a decay interrupted by a 
rather sharp rebrightening, followed by a final steepening. 
We choose to fit the light curve with a Beuermann function 
(Eq.[l| [Beuermann et al.|1999 1 for the first part down to the 
rebrightening: 

KQdll/« 



Fit) = A/[{t/t^ 



{t/tt 



(1) 



where A is a normalisation constant, ar(d) is the slope of 
the rise (decay) phase, and k is a smoothness parameter. 
The time at which the curve reaches its maximum is tmax ~ 
tb(-ar/ad)'/[''^"''-"'^l. 

We model the rebrigthening at about 0.6 d with a simple 
step function. We also tried more complex models by using a 
pulse following the functional form discussed in |Norris et al.| 
( |1996[ ) and in |Guidorzi et al.| ( |2007a| . The quality of data did 
not allow us to study the transition in detail; we report the 
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Table 1. Optical/NIR observations of GRB071010A. The reference time toRB is Oct. 10, 2007 at 03:41:12 l |Moretti et al.|2007| . Data 
are not corrected for dust absorption. For the TAROT data, the "clear" magnitudes are calibrated against the R filter. Data are sorted 
according to wavelength and according to time for each filter. 



Mean date 
(UT) 
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day) 


Exp time 
(s) 
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nass 
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Instrument 


M 


ag 


lit 
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2007 
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11 
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88240 
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1 


2 


U 


VLT + FORS2 


20 


77 


± 





07 


2007 


Oct 


12 


01151 


1 


85790 
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1 


1 


U 
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22 


11 


± 





13 


2007 
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11 


03945 
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60x2 


1 


2 


B 


VLT + FORS2 


21 


16 


± 
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2007 
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12 
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1 
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60X5 


1 


1 


B 
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22 


75 


± 





05 


2007 
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10 
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3 


9 


V 


NOT + PolCor 


20 


20 


± 





09 


2007 


Oct 


11 
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40 X 1 


1 


2 


V 


VLT + FORSl 


20 


35 


± 





04 


2007 


Oct 


11 
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225 X 1 


1 


2 


V 


VLT + FORSl 


20 
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± 
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2007 


Oct 


11 
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1 


2 


V 
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20 


35 


± 
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2007 


Oct 


11 


04188 
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2 


V 


VLT + FORSl 


20 


36 


± 





04 


2007 


Oct 


11 


04496 





89135 


225 X 1 


1 


2 


V 


VLT + FORSl 


20 


36 


± 





04 


2007 


Oct 


11 


04804 





89443 


225 X 1 


1 


2 


V 


VLT + FORSl 


20 


37 


± 





04 


2007 


Oct 


11 


05135 





89774 


225 X 1 


1 


2 


V 


VLT + FORSl 


20 


39 


± 





04 


2007 


Oct 


11 


05441 





90080 


225 X 1 


1 


S 


V 


VLT + FORSl 


20 


40 


± 





04 


2007 


Oct 


11 


05748 





90387 


225 X 1 


1 


S 


V 


VLT + FORSl 


20 


40 


± 





04 


2007 


Oct 


11 


06056 





90695 


225 X 1 


1 


2 


V 


VLT + FORSl 


20 


42 


± 





04 


2007 


Oct 


11 


20925 


1 


05564 


180 X 1 


1 


7 


V 


Keck_I + LRIS 


20 


62 


± 





04 


2007 


Oct 


11 


21238 


1 


05877 


180 X 1 


1 


7 


V 


Keck_I + LRIS 


20 


62 


± 





04 


2007 


Oct 


12 
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1 


87321 


60x2 


1 


1 


V 


VLT + FORS2 


21 


94 


± 





04 


2007 


Oct 


13 


00848 


2 


85487 


60X2 


1 


1 


V 


VLT + FORS2 


22 


75 


± 





04 


2007 


Oct 


10 


20107 





04745 


100 X 1 


1 


7 


R 


Keck_I + LRIS 


18 


15 


± 





06 


2007 


Oct 


10 


20351 





04990 


100 X 1 


1 


7 


R 


Keck_I + LRIS 


18 


25 


± 





04 


2007 


Oct 


10 


28570 





13209 


120 X 1 


2 


2 


R 


Keck_I+LRIS 


18 


94 


± 





04 


2007 


Oct 


10 


28884 





13523 


30 X 1 


2 


S 


R 


Keck_I + LRIS 


18 


98 


± 





04 


2007 


Oct 


10 


29230 





13869 


SOX 1 


2 


S 
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Keck_I + LRIS 


18 


97 


± 





04 


2007 


Oct 


10 


29375 





14014 


SOX 1 


2 


4 


R 


Keck_I + LRIS 


18 


98 


± 





04 


2007 


Oct 


10 


29530 





14169 


SOX 1 


2 


4 


R 


Keck_I + LRIS 


19 


00 


± 





04 


2007 


Oct 


10 
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2 


S 


R 
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20 


00 


± 





20 


2007 


Oct 


11 


04311 





88950 


60 X 1 


1 


2 


R 


VLT + FORS2 


19 


78 


± 





04 


2007 


Oct 


11 


06735 





91374 


180 X 1 


1 


S 


R 


VLT + FORSl 


19 


S3 


± 





04 


2007 


Oct 


11 


19998 


1 


04637 


SOX 1 


1 


7 


R 


Keck_I + LRIS 


20 


05 


± 





04 


2007 


Oct 


11 


20162 


1 


04801 


SOX 1 


1 


7 


R 


Keck_I + LRIS 


20 


01 


± 





05 


2007 


Oct 


11 


20347 


1 


04986 


45 X 1 


1 


7 


R 


Keck_I + LRIS 


20 


03 


± 





04 


2007 


Oct 


11 


20582 


1 


05221 


45 X 1 


1 


7 


R 


Keck_I + LRIS 


20 


04 


± 





04 


2007 


Oct 


12 


02961 


1 


87600 


60 X 1 


1 


1 


R 


VLT + FORS2 


21 


34 


± 





04 


2007 


Oct 


13 


01453 


2 


86092 


60X5 


1 


1 


R 


VLT + FORS2 


22 


14 


± 





04 


2007 


Oct 


14 


00919 


3 


85558 


60x5 


1 


1 


R 


VLT + FORS2 


22 


87 


± 





05 


2007 


Oct 


16 


01841 


5 


86480 


60X 15 


1 


2 


R 


VLT + FORS2 


23 


80 


± 





08 


2007 


Oct 


16 


20948 


6 


05586 


Sx 150 


1 


7 


R 


Keck_I + LRIS 


24 


03 


± 





05 


2007 


Oct 


10 


82711 





67350 


0.1X6000 


3 


2 


I 


NOT+PolCor 


18 


81 


± 





05 


2007 


Oct 


11 


04457 





89096 


60 X 1 


1 


2 


I 


VLT + FORS2 


19 


14 


± 





03 


2007 
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12 


03223 


1 


87862 


60X2 


1 


2 


I 


VLT + FORS2 


20 


65 


± 





03 


2007 


Oct 


14 


01544 


3 
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60x5 


1 


1 


I 


VLT + FORS2 


22 


23 


± 





05 


2007 


Oct 


16 


01565 


5 
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60X 15 


1 


1 


I 


VLT + FORS2 


23 


06 


± 





10 


2007 


Oct 


27 
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le 
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1 


2 


I 


VLT + FORS2 


24 


41 


± 





29 


2007 
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11 


01443 


3C 
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180X9 


1 


5 


I 


VLT + FORS2 


24 


44 


± 





25 


2007 


Oct 


10 
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00381 


10X5 


2 


S 


J 


REM + REMIR 


14 


72 


± 





18 


2007 


Oct 


10 
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01137 


SOX 5 


2 


5 


J 


REM + REMIR 


15 


08 


± 





13 


2007 


Oct 


10 


16737 





01376 


240 


2 


5 


J 


GROND 


15 


16 


± 





11 


2007 


Oct 


10 


17856 





02495 


60X5 


3 





J 


REM + REMIR 


15 


64 


± 





13 


2007 


Oct 


10 


26426 





11053 


30X9 


2 





J 


GEMINI_N + NIRI 


16 


90 


± 





05 


2007 


Oct 


11 


15218 





99856 


1200 


2 


2 


J 


GROND 


17 


76 


± 





10 


2007 


Oct 


13 


01086 


2 


85726 


60X 15 


1 


1 


J 


NTT + SofI 


19 


90 


± 





08 


2007 


Oct 


10 


15519 





00158 


10X5 


2 


S 


H 


REM + REMIR 


14 


43 


± 





15 


2007 


Oct 


10 


16038 





00677 


SOX 5 


2 


4 


H 


REM + REMIR 


14 


02 


± 





10 


2007 


Oct 


10 


16737 





01376 


240 


2 


5 


H 


GROND 


14 


38 


± 





05 


2007 


Oct 


10 


17049 





01688 


60X5 


2 


8 


H 


REM + REMIR 


14 


52 


± 





10 


2007 
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10 
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60X 12 


2 


1 


H 


GEMINI_N + NIRI 


16 


12 


± 





05 


2007 


Oct 


10 


82059 





66698 


60X30 


2 


1 


H 


TNG + NICS 


16 


62 


± 





04 


2007 


Oct 


11 


15218 





99856 


1200 


2 


2 


H 


GROND 


16 


94 


± 





05 


2007 


Oct 


13 


04876 






60 X 20 








NTT + SofI 
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07 


2007 


Oct 


10 


15632 





00271 


10X5 
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S 


K 


REM+REMIR 


13 


42 


± 
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2007 


Oct 


10 
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REM + REMIR 


13 


39 
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10 
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240 
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results requiring the minimum number of free parameters. 
A general discussion is given in Sec.|4] 

The later-time evolution of the afterglow shows a final 
steepening, and this was also modelled with a Beuermann 
function (Eq.jl]). Both for the onset and the final steepen- 
ing the transition from the first to the second power-law be- 
haviour was remarkably sharp and the Beuermann's smooth- 
ness parameter had a high value (k 10, frozen in the fits). 



making the model essentially consistent with what would be 
obtained by using simple power-law segments. The tempo- 
ral decay shown by the Swii/t-XRT data was modelled with 
the same functional form as the optical/NIR where quasi- 
simultaneous data were available. 

The data were fitted both in the spectral and tempo- 
ral domains. The optical/NIR afterglow spectrum was mod- 
elled with a simple power law. Galactic dust absorption was 
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Table 2. Best-fitting parameters. 



Figure 1. The field of GRB 071010A at about 3.8 d after the GRB 
in the R band observed with the VLT equipped with FORS 2. A 
nearby source at ij = 23.22 it 0.06 at about one arcsec Eastward 
of the afterglow is visible (see Sect.jsjl. The seeing of the image is 
~ 0.7 arcsec. 



removed and rest-frame absorption was modelled with the 
absorption curves of the Milky Way (MW), the Large and 
Small Magellanic Clouds (LMC, SMC) (|Pei||1992[ ), and for 
a starburst galaxy (SB) (Calzetti et al. 20001. For the I 



band we also added a constant component, possibly the host 
galaxy or a supernova (SN), as late-time data (later than 
about 10 d) show a clear flattening. However, the presence 
of a red point-like object ~ 1" east of the afterglow, with 
R = 23.22 ± 0.06 mag and / = 22.31 ± 0.07 mag, makes 
late-time photometry less reliable (see Fig.[T]). This source 
could also contribute to the NIR flux at about 2.8 d, when 
the observing conditions did not allow a reliable separation. 
A simple extrapolation to the NIR bands for R — Kml mag 
would suggest that it could be comparable in brightness to 
the afterglow. We therefore did not consider these data in 
the fits. 

The model gives an acceptable fit (Xdof — 114.1/90 ~ 
1.27), particularly if one considers the number of telescopes 
contributing the dataset under different observing condi- 
tions, filters, etc. A summary of the best-fitting parameters 
is reported in Table[2] 

The early-time observations by TAROT, REM, and 
GROND allowed us to identify a rising phase, a maxi- 
mum, and then a regular decay phase. The fit yields tn 



420^ 



s (~ 7min) 



-0.88 



-1-0.43 



and Qdi = 0.71 



+0.03 



for the time of the maximum and the rising and decay- 
ing power-law indices, respectively, in substantial agreement 



with the findings reported by Klotz et al. ( j2007b ). The data 
do not suggest the presence of any spectral evolution. 

A rebrightening peaking at about 0.6 d dominates the 
afterglow at intermediate times. The required flux ratio in- 
crease (~ 1.9) is remarkably high and, within the uncer- 
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0.2i: 
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tainties, the same in the optical/NIR and X-rays. The later 
evolution of the afterglow requires a break, consistent with 
being fully achromatic in the optical/NIR and X-ray bands, 
at ibroak = 0.96lQ Qg d, with a post-break decay index of 
Qd2 ~ 2.071q q7, followed by a final flattening likely due to 
the contribution from the GRB host galaxy. The best-fltting 
model is shown superposed on the data in Fig. [2] 

The entire afterglow spectral evolution is consistent 
with being achromatic. In the optical/NIR the broad-band 
spectral energy distribution index is (3 = 0.76!to'26 
sizeable rest-frame extinction Eb-v = 0.21 
having removed the Milky Way absorption along the line of 
sight. Among the various extinction curves, only that of the 
SMC gave satisfactory results. The r edshift of GRB 071 010A 
would move the 2175 A bump (see Fynbo et al.|2007 for a 
discussion about observations of this feature in GRB after- 
spectra), typical of Milky Way absorption ( Pei||1992 |, 



with a 
ao5 mag after 



into the B band, therefore easily detectable, if present, in 
our data sets. 

The X-ray spectrum was modelled assuming a simple 
power law with neutral absorption in the Milky Way and 
in the host galaxy. Analysis of the hardness ratio showed 
that there are no significant spectral variations during the 
X-ray observations. The spectral analysis of the X-ray data 
alone (Xdof — 14.4/19 « 0.76) provided satisfacto ry results 



Fixing the Galactic absorption at 6.5 x 10 cm ( Dickey & 



Lockman 1990) the rest-frame absorption and the spectral 
index turned out to be 1.74(+o 37) x 10^^ cm"^ and f3 = 
l-46to ig, respectively. The optical/NIR and X-ray spectral 
indices are consistent within the uncertainties with /3x ~ 
/3opt+0.5, i.e. with the cooling-break between the two bands. 
The broad-band SED is shown in Fig.|3] 

The prompt time-averaged spectrum of GRB 071010A 
in the Swift BAT 15-150 keV energy band is consistent with 
a single power law with photon index significantly higher 
than 2 ( [Krimm et al.|[2007[ ), suggesting that the peak en- 
ergy, -Bpeak, of this cveut is close to, or below, 20-30 keV. In 
order to better investigate this issue, we reduced and ana- 
lyzed the BAT time-averaged spectrum. The fit with a sim- 
ple power law yelds a photon index of 3.1±0.7, conflrming 
that GRB 071010A is a soft event, an X-ra y flash (XRF). By 
fitting the spectrum with a Band function | |Band et al.|1993"| 
with a fixed at —1, we derive an upper limit to the peak en- 
ergy iJpcak of ~35 keV, which corresponds to ~70 keV in the 
GRB cosmological rest frame. Then, by varying Bpeak ftom 
its upper limit to and adopting the method of |Ghirlanda| 
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et al. (20041 and Amati (20061, we estimated an isotropic- 



equivalent radiated energy in the 1-10000 keV cosmological 
rest-frame energy band, -Eiso, of (3.6 ± 2.9) x lO^'^ erg. 



4 DISCUSSION 
4.1 The onset 

The temporal resolution of the observations does not allow 
us to check in detail the early-time density profile of the cir- 
cumburst medium or possible chromaticity in the evolution. 
The modest duration of the prompt event and the late onset 
time argues against a scenario with frequent multiple peaks 
as for GRB 070311 ([Guidorzi eFal. 2007a). A smooth rising 



phase can also be produced by a decreasing extinction in the 



case of a radially decreasing circumburst density ( Rykoff et 
al. 2004[ ). However, the detection of the afterglow in the 



U band (Table[TJ eliminates this last possibility due to the 
high absorption required. An alternative giving a chromatic 
maximum could be the passage of the typical synchrotron 
frequency in the optical/NIR bands. The time spread for the 
maximum would depend on the width of the observed fre- 
quency band, ti/t2 oc {iyi/i'2)~^^^ , with the higher frequen- 
cies peaking first. The passage would also produce a spectral 
change from a positive (1/3) to negative [— (p— 1)/2] spectral 
index in the slow-cooling case ( |Sari et al.|1998[ [Chevalier fc 
pl|p99l ). A gain, this is not supported by the data. 

The maximum during the early afterglow evolution can 
then be interpreted as the afterglow onset as for the case of 
GR B 0604 18 and GRB060607A ( [Molinari et al.|2007|[jin¥ 
Fan 20071. Following the formalism reported in Eqs. 1 and 



2 of Molinari et al. (20071 in the thin-shell case since the 



prompt duration was much shorter than the peak of the 
afterglow emssison, we can derive, for GRB 071010A, Fo ~ 
150 for a uniform interstellar medium (ISM) and Fo ~ 40 in 
the case of a wind-shaped density profile. These are lower 
than the values inferred for GRB 060418 and GRB 060607A 
due to the late occurrence of the maximum, but still within 
the theoretical expectations for the external-shock scenario 



( Sari & Piran 1999 Lithwick & Sari 2001 Zhang et al. 2006 1 



The observed maximum cannot be related to reverse 
shock emission. The reverse shock should be in general 
a short-lasting phenomenon, and in the thin-shell case it 
should peak slightly before the deceleration time, roughly 
just before the afterglow onset, and dominate it in case the 
typical synchrotron emission is close enough to the optical 
band ( [Kobayashi fc Zhang|20"07| | Jin fc Fan|2007|). The rapid 



decay (/ ~ t~''''^, Sari fc Piran| 1999 Kobayashi fc Zhang 
2007[ ) predicted for the reverse shock emission is inconsistent 
with the much milder decay observed after the optical flux 
peak (Table[2| 



4.2 The rebrightening 

The occurrence of rebrightenings or flares during the tem- 
poral evolution of GRB afterglows is not uncommon. The 
long time coverage usually provided by the Swift XRT af- 
fords the study of a large sample of flares from the spectral 

or temporal and mor- 



and energetic (Falcone et al, 
phological ( [Chincarini et al 



20071 



20071 points of view. In the 



common, while in a remarkable fraction of cases a rebright- 
ening due to discrete energy injection was included in the 
modelling ( Johannesson et al.|2006 |. 

Our dataset does not permit a detailed analysis of the 
temporal and spectral evolution of the rebrightening around 
0.6 d. The simplest approach is to model the transition as 
due to the discrete injection of a sizeable amount of en- 
ergy in the flreball interacting with the external medium by 
means of refreshed shocks. The flux density produced by an 
external shock scales with the energy content of the fireball 
apart from minor differences between the optical/NIR and 
X-ray bands ( Panaitescu fc Kumar|2d00 l which are not de- 
tectable in our data. Therefore, to account for the observed 
/inj// ~ 1.9, an amount of energy comparable to the energy 
content of the fireball needs to be supplied to the system. 
Remarkably, a discrete episode of energy injection is not 
supposed to modify the spectrum of the synchrotron radia- 
tion emitted at the shock front if the cooling frequency does 



not enter or cross the observing bands ( Panaitescu fc Kumar 



2000 1 , in agreement with the achromatic evolution of the af- 



terglow of GRB 071010A. A discrete energy injection, com- 
parable to the energy content of the flreball, was also sug- 



gested by de Ugarte Postigo et al. (20071 for GRB 050408 



Also in the case of the XRF 050824 a rebrightening was ob- 
served at rather late time and a large delayed energy in- 
jection was invoked ( SoUerman et al!||2007 l. In principle, a 
rebrightening could be induced by a density variation in the 



circumburst medium as proposed by Lazzati et al. ( 2002 1 for 
the GRB 021004 (but see also |Nakar fc Granot|2007| ). How- 
ever, this explanation cannot be applied to our case, since 
the X-ray range is above the cooling frequency and there- 
fore insensitive to density variations (see, e.g. 



Sari et al. 



19981. A transition from wind to homogenous medium, due 



to the crossing of the wind reverse shock could also generate 
a rebrightening, as discussed in Pe'er fc Wijers ( |2006[). Other 
explanations are also possible, such as a double jet ( Bcrger et 
al.|2003 l or a patchy shell ( |Rees fc MeszaroS| 1998,|Meszaros 



et al. [1998 1. A structured jet with a bright narrow core and 
fainter wings can also produce a rebrightening close to the 



jet-break time (Salmonson 20031, though in such a case it 



would only be due to perspective rather than to real energy 
injection, possibly driving to a sharp transition if seen from 
a suitable angle. 

Finally, although not required by the data, the temporal 
proflle of the rebrighteiung is also consistent with that of a 
fast rise and exponential decay pulse simply superposed on 
the underlying afterglow. The relatively long duration of the 
pulse, as compared to the time of occurrence, At/t ~ 1, is 
compatible with an origin of the emission at the external- 
shock radius (see discussion in Guetta et al.|2007 Guidorzi 
|et al.|2007al ). 



4.3 Late-time afterglow, light-curve break, and 
SED analysis 

Following the rebrightening at about 0.6 d, the afterglow 
evolution shows a sharp break at about 1 d consistent 
with being achromatic in the X-ray and optical/NIR bands 
(Fig.[2]|. FoUowing |Liang et al.| (|2007| only seven GRB af- 



optical/NIR the occurrence of large flares seems to be less 



terglows have been found so far with such an achromatic 
transition. Moreover, the general temporal behaviour of the 
optical/NIR and X-ray light curves seems to be the same, 
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Figure 2. The optical, NIR, and X-ray light curves of the afterglow of GRB 071010A. The dashed lines show the fit to the temporal 
decay of the early and late-time afterglow. Fit parameters are reported in Sect.js] On the left we show the whole light curve fitted with 
the adopted model. On the right the epoch around Id is magnified. At late times, the long-dashed line shows the possible contribution 
to the / band by a supernova component using SN 1998bw i jGalama et al.|[l998[ l as a template, at the redshift of GR,B071010A. The 
effect of the absorption in the optical is included. 



supporting at least qualitatively the possibility that we are 
observing emission from a single outflow with a transition 
due to a beamed geometry. 

The standard external-shock model predicts specific re- 
lations between temporal and spectral parameters. However, 
it is a well-established observational fact that GRB after- 
glows show more complex behaviours than expected before 
the launch of Swift. Many mechanisms have been proposed 
to interpret the new features discovered with Swift: de- 
layed energy injection, late central engine activity, refreshed 
shocks, high-latitude emission, etc. (e.g., [Zhang fc Meszaros| 
2004[ [Zhang et al.|2006| ). Most of these ingredients affect the 
light curve shape in a number of ways. However, since the 
fundamental emission process in the soft X-rays and in the 
optical/NIR is still synchrotron, the spectra should follow 
the general predictions of the external-shock scenario. We 
therefore begin the analysis starting from the broad-band 
SED from the optical/NIR to the X-rays. The evolution of 
the afterglow, including the early phases, is consistent with 
being fully achromatic. We derive the SED at about 0.8 d in 
order to minimise the need to extrapolate the observations 
adopting a temporal behaviour (Fig.|3|. 

The amount of rest-frame neutral absorption required 
by the X-ray data {Nh « 1.7 x lO^^cm"^) is within the 
range observed for a set of well-studied Swtft GRB X-ray 
afterglows ( [Campana et al.]|2006[ ). In the optical/NIR the 
observed SED curvature clearly indicates the need for a sub- 
stantial absorption, as derived in Sect. [3] with Eb-v ~ 0.21 
mag. Such an amount of intrinsic absorption in the opti- 
cal/NIR is uncommon among GRBs with sufficient data 



quality to allow a meaningful study (see Stratta et al. 2005 



[Tagliaferri et al.|2006||Starhng et al.|2007||Kann et al.|2007 

for a comprehensive analysis). As for most of the GRB opti 
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Figure 3. SED of the afterglow of GRB071010A about 0.8 d 
after the burst. The dashed line is a fit to the observed data in 
the optical/NIR. The solid line refers to the intrinsic SED once the 
effect of dust absorption in the Milky Way and in the GRB host 
galaxy is removed. The X-ray and optical/NIR broad-band SED is 
consistent with the predictions of the synchrotron external-shock 
scenario with a cooling break between the optical/NIR and the 
X-ray and /3x = /3opt + 1/2. 
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cal afterglows, a SMC-like extinction curve seems to provide 
a better fit to the data ( |Kann et aL]|2006[ [2007| , although 
the gas-to-dust ratio Nh/ Eb-v is higher than in the SMC, 
again in agreement with the findings of Starling et al. ( 2007 1. 



The spectral indices in the optical/NIR and X-ray 
bands allows us to constrain the electron energy distribution 
index p {N oc 7"*"). With the cooling frequency in between 
the optical/NIR and the X-rays (probably just below the X- 
ray range, as shown in Fig.[3|, we have /3opt ~ (p— l)/2 and 
/3x =p/2 and therefore popt = 2.52+° g^ andpx = 2.92+° gg, 
independent of the circumburst matter-density profile in the 



slow-cooling regime ( Zhang fc Meszaros^2004^ and references 
therein). These values for p, in particular from the X-rays, 
lit in the high tail of the distribution derived by |Tagliaferri| 
et al. (20061 or Shen et al. (2006), although still consistent 



with it. The value derived from the optical is also consistent 
with the typical results {p « 2.2) of numerical simulations 
of relativistic shocks ( Achterberg et aL]|2001 Vietri 20031. 
Knowing the electron distribution index we can predict the 
decay rate in the optical/NIR before the occurrence of a 
jet break, considering also the constraints from the X-rays: 
Qpre.opt « 1.2 - 1.5 (0.75-1.5 from the optical/NIR only) 
in the case of a constant-density radial profile of the cir- 
cumburst matter, or Qpro.opt ~ 1.7 — 2.0 (1.25-2.0 from the 
optical/NIR only) in the case of an environment shaped by a 
wind from a massive star. For the X-rays (frequencies above 
the cooling frequency) we have apre,x ~ 1.4 — 1.75 (1.4-2.0 
from the X-ray only) independent of the density profile (and 
unconstrained by our data). 

It is clear from the fit results, reported in Sect. [3] and 
shown in Fig.[2j that interpreting the break at ~ 1 d as a jet 
break, the optical/NIR decay before the transition would be 
in essential agreement with the expectations for a constant- 
density medium only for the p value derived from the optical. 
Otherwise, a continuous additional injection of energy as re- 
quired by the modeling of numerous Swift X-ray afterglows 
( Panaitescu et aL||2006 1 should be considered. After the jet 
break, the decay index in the slow-cooling regime must be 
the same at all frequencies (higher than the typical syn- 
chrotron frequency, as expected ~ 1 d after the high-energy 
event), and independent of the circumburst matter-density 
profile. For p > 2, as in our case, it is simply cipost = P- Again 
the observed decay indices are consistent with the expecta- 
tions only with the value derived from the optical/NIR. 

In principle, the post-break decay indices would also 
be consistent with a wind model before the occurrence of 
a jet break with p ~ 3.1. This p value would be better 
supported by the X-ray spectral data, although the achro- 
maticity of the transition at ~ 1 d is not consistent with this 
scenario. The distribution of p values, although likely incon- 
sistent with a unique p ~ 2.2 ( Tagliaferri et aL]|2006 Shen 
et al.|2006 l, supports p < 3 for essentially all the afterglows 
once the analysis is based on the spectral parameters rather 
than on the temporal behaviour. 

Summing up, the hypothesis that the sharp achromatic 
transition at about 1 d be a real jet break appears consis- 
tent with the data if we assume that continuous energy 
injection is affecting the temporal decays. This would be 
one of the very few instances of a jet break observed both 
in the optical/NIR and the X-rays. The consistency of the 
GRB 071010A parameters with the Ghirlanda relation (see 



below) lends some indirect support to the jet-break inter- 
pretation. 

After about ten days the afterglow light curve flat- 
tens considerably, plausibly as a result of the emergence of 
the host-galaxy light. Our data quality does not permit to 
draw any conclusion about the possible presence of a SN 
component affecting the light curve and mimic the host- 
galaxy contribution (see Fig.[2|). A SN component compa- 
rable to SN 1998bw ( |Galama et"aL]|1998[ ) at the redshift of 
GRB 071010A and including the effect of the fitted absorp- 
tion in the optical would be in the I band about 0.8 mag 
fainter than our last measured point. 



4.4 Amati and Ghirlanda relations 

It is well known that the position of a GRB in the radiated 
energy vs. spectral peak photon energy plane can provide 
useful clues on the nature, physical origin, and geometry 
of its emission. Indeed, while all "normal" long GRBs and 
XRFs show a clear correlation between these two quanti- 
ties, short GRBs and peculiar sub-energetic GRBs (980425 
and, possibly, 031203) show a different behavior (e.g.. Am- 
ati|2006 1. In addition, coUimated GRBs seen off-axis are ex- 



pected to deviate from these spectrum-energy correlations. 

The -Epeak and -Eiso values derived from the analysis 
of the time-averaged spectrum (see Sect.jsJ are fully con- 
sistent with the -Bp,i — -Biso (Amati) relation ( [Amati et al. 
|2002 Amati 2006 1 , showing that the prompt emission of 
GRB 071010A is not peculiar and indicating that, if the 
emission is collimated, the off-axis angle is small. In addi- 
tion, with a break interpreted as a jet break at about 1 d and 



by assuming, following Ghirlanda et al. (20041, a constant 
circumburst matter-density profile (n = 3 cm~'^) and a fire- 
ball kinetic to radiated energy conversion efficiency of 0.2, we 
derive a jet opening angle of 8.2±1.1 deg and a coUimation- 
corrected energy of (3±2) x 10"*^ erg. Thus, GRB 071010A is 
also fully consistent with the i?p,i — (Ghirlanda) relation 
( Ghirlanda et al.,2004}|Nava et al.|2006| . This evidence fur- 
ther supports that the possibly achromatic break observed 
in the afterglow light curve is due to collimated emission. 

Finally, if we fit the time-averaged spectrum of 
GRB071010A with a cut-off power law with index fixed at 
— 1, we can constrain the cosmological rest-frame peak en- 
ergy to 32+21 keV. This value, combined with the values of 
_Eiso and E-y derived above, make the representative point 
for GRB 071010A consistent with the best-fit power law of 
both the Ep^i — Eiso and Ep,i — E~, relations, thus reinforcing 
the above considerations. 



4.5 Polarisation 

Almost coincident with the achromatic break, we obtained 
polarimetric observations, providing a rather robust upper 
limit (0.9% at 95%, 1.3% at 3(t) suggesting null polarisa- 
tion. Polarimetric measurements of GRB optical afterglows 
have been performed for several events ( Covino et al.|2004 |. 
In general, the polarisation level was moderate (2-3%), al- 
though in a few cases evidence for variations related to the 



afterglow evolution was found (Lazzati et al 



^t al^i2003 



2003 



Bersier 



Rol et al. 112003 


IGreiner et al.||2003| |Gorosabel 


|Lazzati et al. 


2004 


1. The simple detection of 
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variable polarization, intrinsically related to the afterglow 
emission, is still one of the most important observational 
tests in favour of the external relativistic shock model with 
physical beaming, where some degree of polarisation is nat- 



urally expected (Malesani et al. 
Lazzati|2006l ) 



2005 Covino et al. 2005 



Almost all of the observations have been performed with 
the goal of disentangling the geometry of the outflow using 
time-resolved polarimetry. After an initial period of enthusi- 
asm, the required intense and delicate observational efforts 
and the high degree of complexity of the afterglow light 
curves provided by Swift cast some doubts on the actual 
possibility of deriving meaningful constraints from simple 
models ( [Ghisellini fc Lazzati,1999||Sari|1999l ). These proved 
to be applicable only to smooth, regular, afterglow light 
curves. Nevertheless, Rossi et al. (20041 showed that the 



polarisation degree and position-angle evolution is strongly 
dependent on the assumed jet structure. At least in a cou- 
ple of cases a homogeneous jet geometry could be formally 
ruled out: i.e. for GRB 020813 ([Lazzati et al.][2004| and for 
GRB 030328 ( [Maiorano et al.|2006[ ). 

The diagnostic power of time-resolved polarimetric ob- 
servations relies on the strong coupling between the assumed 
jet geometry and polarization signature. For homogeneous 
jets there should be two polarization maxima, before and 
after the jet-break time, with a rotation of the position an- 
gle by 90°. For structured jets, (i.e., jets with a more ener- 
getic core and fainter wings), the maximum of polarization 
should almost be coincident with the jet-break time, and the 
polarization angle constant. Therefore, a polarimetric mea- 
surement almost coincident with a jet break, as in the case 
for GRB071010A, is of remarkable value. The very low up- 
per limit that we have obtained (Sect.jsjl argues against a 
structured jet. However, the lack of time-resolved polarime- 
try through the evolution of the afterglow prevents us from 
drawing firm conclusions. The polarised fraction for the ra- 
diation emitted by an afterglow depends basically on the 
degree of asymmetry of the system, which in turn depends 
on the line of sight of the observer with respect to the jet 
symmetry axis. With a single observation the possibility of a 
very small offset angle cannot be ruled out. Setting aside the 
jet-break hypothesis, it is still true that, whichever physical 
process is responsible for the (achromatic) steepening ob- 
served in the afterglow at about 1 d, it generated radiation 
with only a very low level of polarization. The same is true if 
we model the rebrightening at about 0.6 d with a fast rise ex- 
ponential decline pulse dominating the afterglow luminosity 
for the pulse duration. 



5 CONCLUSION 

Our campaign on GRB 071010A is an example of extensive 
monitoring of a GRB afterglow in the optical/NIR. This 
allowed us to detect an initial rising phase which can be 
interpreted as the afterglow onset. The peak time is rather 
late, about 7min after the GRB. This would directly imply 
a relatively small initial Lorentz factor although still within 
the fireball model expectations. Interesting features are the 
rebrightening at ~ 0.6 d followed by an achromatic break 
~ 1 d after the burst. The steepening is observed both in 



the optical/NIR and in the X-rays bands, a feature seldom 
observed even in the Swift era. 

The spectral and temporal relations before and after 
the break are roughly consistent with the predictions of the 
afterglow synchrotron external shock model with a cooling- 
break in between the optical/NIR and X-ray bands and as- 
suming also continuous energy injection. The rebrightening 
episode can be interpreted as being due to the occurrence of 
a refreshed shock injecting energy in the system comparably 
to the original energy content of the fireball. Other expla- 
nations are still possible, relating the rebrightening to the 
density profile of the circumburst medium or to structured 
jets. The optical/NIR spectral energy distribution requires 
absorption in the host-galaxy satisfactorily modeled with 
the SMC extinction curve. 

GRB 071010A is consistent with the Amati relation and, 
if the achromatic steepening is interpreted as a jet break, it 
is consistent with the Ghirlanda relation too. Independent 
of the emission process, the entire afterglow evolution is con- 
sistent with being independent of color from about 2 min up 
to a few days. The SED is consistent with synchrotron emis- 
sion and requires rest-frame absorption both in the X-rays 
and in the optical/NIR. 

We also performed a polarimetric observation almost in 
coincidence with the occurrence of the achromatic break. A 
tight (0.9%) upper limit was derived. In the absence of time- 
resolved polarimetry, we cannot rule out that this afterglow 
was characterized by a low level of polarization throughout 
its evolution. 

Finally, the rich dataset of this burst and the simul- 
taneous availability of multi-wavelength data is a positive 
example of a new approach we are trying to apply to GRBs 
localized by Swift, favouring more complete coverage of a few 
well-studied GRBs. In addition, the observation of an initial 
rising phase, a rebrightening episode, and an achromatic 
break, as well as the availability of a well-sampled broad- 
band energy distribution, could offer a valuable workbench 
for alternative interpretative scenarios such as the "cannon- 
ball" pe Rujula||2007[ [Dado et alT]|2007[), the "fircshell" 



( Bianco et al.|2007[ ), and the "quark nova" ^Staff et al. 2007| 
models in order to reach a better understanding of GRB af- 
terglow phenomenology. 
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